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COMPRESSIBLE POTENTIAL FLOW WITH CIRCULATION ABOUT A CIRCULAR CYLINDER 

By Max. A. Hkabijst 



SUMMARY 

The 'potential function for flow, with circulation, of a com- 
pressible fluid about a circular cylinder is obtained in series 
■form including terms of the orders of M* where Mis the Mach 
number of the free stream. The resulting equations are used to 
obtain pressure coefficient as a function of Mach number at a 
point on the surface of the cylinder for different values of circula- 
tion. The coefficients derived are compared with the Qlauert- 
Prandtl and Kdrmdn-Tsien approximations which are func- 
tions of the pressure coefficients of an incompressible fluid. 
For the cases considered, the values of the pressure coefficients 
computed from the theory were found to be somewhere between 
the two approximations, the first underestimating and the 
second overestimating it. 

INTRODUCTION 

In the two-dimensional irrotational flow of a compressible 
fluid, where the expansion is assumed to be adiabatic, the 
velocity potential is known, to satisfy a nonlinear partial 
differential equation of the second order. For subsonic 
velocities, at least three methods are known for the approxi- 
mate solution of this equation. They are usually denoted 
as the method of small perturbations, the Rayleigh-Janzen 
method, and the hodograph method. 

The method of small perturbations (references 1 and 2) 
assumes that velocity changes which are brought about by 
the airfoil in the uniform parallel air stream are small in 
comparison with the velocity of the undisturbed stream. 
Under this assumption it is possible to introduce new vari- 
ables which reduce the differential equation to a Laplace 
equation, and, as a consequence, the problem becomes one 
concerning flow in an incompressible fluid, provided the body 
is assumed distorted to correspond to the change of variables. 
The assumed distortion consists in expansion of the dimen- 
sions of the airfoil perpendicular to the direction of the free 
stream in the ratio l/^/l-M 2 , where Mis the Mach number 
of the undisturbed stream. 

The Rayleigh-Janzen method (references 3 and 4) assumes 
that the general expression for velocity potential may be 
written as a series in rising powers of M and with variable 
coefficients. These coefficients can be shown to satisfy 
certain Poisson differential equations and, if the equations 
are integrable, the solution becomes a matter of determining 
these coefficients. Successive steps, however, become in- 



creasingly laborious and the convergence of the series may 
be slow, even at relatively small Mach numbers, if the shape 
of the body is such that the speed of sound is approached 
locally. Solutions, using this method of attack, have been 
carried out by C. Kaplan (references 5 and 6), S. G. Hooker 
(reference 7), I. Imai (reference 8), K. Tamada and Y. Saito 
(reference 9), and L. Poggi (reference 10). Poggi introduced 
certain refinements and some of the preceding references 
employ this process. It is tantamount to using the so-called 
Neumann function in solving given Poisson equations and 
will be discussed in the appendix. 

The hodograph method is ascribed by writers on that 
subject to P. Molenbrock and A. Tschaplygin. Instead of 
expressing the velocity potential as a function of coordinates 
in the Cartesian or polar plane, the magnitude of velocity 
V and its inclination 0 to an assumed axis are chosen as 
independent variables. The resulting 'differential equation 
is linear and can be further simplified by replacing the pres- 
sure-volume relationship for adiabatic expansion by the 
equation of a line tangent at a point corresponding to the 
state of the fluid in the ambient stream. This artifice was 
suggested by T. von Karman (references 2 and 11) and used 
successfully by Hsue-Shen Tsien (reference 12) . K. Tamada 
(reference 13) has also applied Tsien's more general results 
on elliptic cylinders to compressible flow past a circular 
cylinder. 

One noteworthy result of the hodograph method has been 
the Karman-Tsien expression for pressure coefficient P in 
terms of Mach number M and Pu-a, the pressure coefficient 
for M=o. This expression may be written 

It always gives, for negative pressure coefficients, a result 
greater in absolute value than the Glauert-Prandtl formula 
which is based on the method of small perturbations, 

P=Pjf -° Vl-M 2 (2) 

and is currently accepted as the more accurate of the two. 

From equations (1) and (2) it is possible to compute the 
critical Mach number M B , the value of M at which the local 
speed of sound is attained, in terms of P M -o- The relations 
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involving M c and P m -q, corresponding respectively to 
formulas (1) and (2), are 



it-o- 



and 



1+71=17? 2 



M-0- 



2yi-M ( 



(3) 
(4) 



The difficulties inherent in the two latter procedures are 
quite as distinctive as their respective approaches to the 
problem. As stated before, the Rayleigh-Janzen method 
employs classical mathematics, the required terms being 
solutions of Poisson equations with given boundary condi- 
tions, but the work involved is arduous. In the hodograph 
method the principal difficulty is to determine proper bound- 
ary conditions in the V, 6 plane. In available calcula- 
tions the solution is given with a slight distortion in the 
given boundary. It is possible to correct this distortion, 
in some cases, so that the final results are not too seriously 
affected. When the flow around the body involves circu- 
lation, however, the change in the boundary is more serious, 
for nonperiodic terms appear and the boundary is no longer 
a closed curve. 'Added circulation does not involve any 
essential variations in the Rayleigh-Janzen method, however, 
and in this report the velocity potential for such compressi- 
ble flow about a circular cylinder has been derived. Since 
no theoretical study has been presented, as far as is known, 
to determine the error in the Karman-Tsien pressure coeffi- 
cient, the results obtained in this report furnish a means of 
approaching this problem. The results of such calculations, 
for various values of circulation, are therefore included. 

ANALYSIS 

Consider a gas obeying the adiabatic law and flowing 
irrotationally in two dimensions. Its equation of motion 
may be written in polar coordinates in the form 



where 



$ velocity potential 
b 5 $ . 1 b$ . 1 



y ratio of specific heats of gas 

c 0 velocity of sound in undisturbed flow 

U velocity of free stream 

M Mach number of free stream I 



V s local velocity squared +p (^)*J 



By the introduction of the variables $ and v so that 

■ * a V 
4> = jj ana p=£j 

equation (5) may be written in the form 
where 

Following the method of Rayleigh and Janzen, assume that 
$ may be developed in a series of ascending powers of M 2 
so that 

4>=<p 0 +M'4> 1 +M*<h+ (8) 

After substitution of equation (8) in equation (7), elementary 
calculations show that 



where 



v i =v*+v 1 1 M i +v» 2 M i + (9) 

In a similar manner, equations (8) and (9) may be substituted 
in equation (6) and on equating coefficients of the same 
powers of M, the following relations for <fo, fy, <fe, . . . 
result: 

VVo=0 (11a) 



s _1 (Ma dfo 2 , 1 dz^A 
V 2 <fe=| (7-D W-l) V 2 &+ 



(Hb) 



1/&&, W , 1 &*o&sA j. 

1 bvp 3 1 dfo bv 0 * \ 
2\br dr^r 1 b9 be J 



(He) 



If the equations (11a), (lib), and (11c) can be solved suc- 
cessively for 4>o, <t>u <fa, • • • the values may be substituted 
in equation (8) to get the potential function for the flow of 
a compressible fluid. A step-by-step procedure is therefore 
established whereby any desired degree of approximation 
to <f> may be obtained, provided the value of M is within the 
region of convergence of the resulting series. Equation (11a) 
is the differential equation satisfied by the potential function 
in the case of ^compressibility. Once this potential func- 
tion is known it is used to evaluate the right-hand member 
of the second equation, the solution of which furnishes the 
second term in the development of <£. The method of ob- 
taining further terms follows the same general procedure. 

Consider now the case of a right circular cylinder of 
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i n fi n ite length in a compressible fluid, the axis of the cylinder 
being at right angles to the direction, of steady flow. In 
determining the velocity distribution about the cylinder, the 
problem may be treated two dimensionally with a circle as 
the boundary curve and the equations established in the 
Rayleigh-Janzen method may be applied directly in the 
following manner. The radius of the circle is arbitrarily 
assumed equal to one, and a polar coordinate system is chosen 
with origin at the center of the circle and polar axis extending 
downstream. The flow about the circle is assumed to be 
that resulting from the combination of uniform stream 
velocity and circulation about the cylinder. Under these 
conditions the classical expression for <f> 0 is well known. 
It may be written 



cos 0- 



where r is the circulation around the circle, measured 
positive in a clockwise direction. For ease of computation it 
is convenient to set 

TrU 



and, as a consequence, 



cos 0—^0 



The boundary conditions, in general, are 



and 



^=0forr=l 



-^=cos 0 for r=«> 



From equations (12) and (10a) 

». a =(l+^)-| cos 20+kQ;+£) sin 0+§J 

This result, together with equation (lib), gives 

2 



(12) 

(13a) 
(13b) 

(14) 



VVi— (tt+p) cos e +p cos 30+ 



cos 0 (15) 



The more elementary methods of integration lead to certain 
difficulties when an attempt is made to solve for <fo, in 
equation (15). These difficulties result from nonperiodic 
terms in the particular integral and resultant trouble in 
determining such constants of integration that the necessary 
periodicity, in terms of 0, is maintained in the final expression 
for the potential function. This difficulty may be obviated, 
however, by established methods. (See appendix.) It 
follows that the solution of 



v 2 a= 



sin m0 



(16) 



satisfying the boundary conditions 



IS 



o sin m$ /(«— 2) ml. , , „ 

Q -m(m-8+2) ^+8-2)1-7= p=s)(^hen m+2^s) 



and 



(17a) 



0 —sin mflT 1.1, 1 , , 

" = mr™ l2m + 2 lo S r )( whei1 m+2=s) (176) 

The veracity of these solutions, together with analogous ones 
existing when sin m0 is replaced- by cos m0, may be checked 
easily by substitution in equation (16). 

Since equation (15) is a linear differential equation, its 
solution is deteimined by considering each term of the right- 
hand member and summing the individual integrals obtained 
by means of equations (17a) and (17b). The final result is 

« 3 * ^ 2 * (-8+6^+24?+^)+ 

*™ (18) 

In the evaluation of <fa the calculation follows the same 
pattern of development. From equation (10b), together 
with equations (12) and (18) 

■^y?-5.+5.)+« » (^-^)+ 

"»(- 4 Vp- 8 ?+5.-^)t ■ 



logr logr 

2r 3 



y^0{ + ^) 



(19) 



This result, together with equations (12), (14), and (18), 
substituted in equation (11c), gives 

V s <fe=(7-l){cos 0 (-^+p-|+^)+cos 30 

cos 50 (=r)+K sin 20 (^-^-^+^) + 

^sin4e(| i 4-| i -^)+X s cos0(=^-|,-|- 7 +| 5 )+ 

i?cos3^4^4j^)+X3sm2^- 8 A +i l 3 ) + 



K* cos 0 (jopt) +{cos 0 (■ 



-32 , 39_15 , 11 \, 



3r* 1 2r 7 



on /19 .3 5 , 1 \ . 
cos 30 {jF+p-zpt+^u ) +cos 59 



\2r> i-sj 



+ 



K sin 20 



(=1 11 
V 4r s_ 3r J 



25 14,. 3 6 logr , 4 logr 



12r« 3r rr 2r 10 



r 8 



) 



+ 
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K sin 40 



(±+11 



, 2 5 1 , 3 log r \ 

i-p I2r 8 " 1 "l2r 10 " t " r* ^ 



™ /-3 2 13 logr 

z» cos e (^--g-^— 



logr\ 
2r 7 y + 



X s cos 35 



_3_ 
"4^ 



JP sin 20 (-J-gp--^ 



1 . log r . 3 log r log r\ 
4r s ~ t ~ 2^ r 5 r 7 y" 1 " 

te)+^cos.(^)l (20) 



To integrate, formulas (17a) and (17b) are again resorted 
to. The method of integration given in the appendix also 
provides integrals corresponding to the new type of terms 
appearing in the right-hand member of equation (20). Thus, 
the solution of 



v%= logr|n_m? 
satisfying boundary conditions 



(21) 



is 



when m^s— 2. 



sin md f 1 . 1 

"27nr» i(m+s-2) ai "(m-«+2) 2 

sin m B f 2m log r 

2mr»- 2 i(m+s— 2) (m— s+2)~ 

1 1 

(m+s-2) 2 (m— a+2) : 

When m=«— 2, 



_ —sin 77i5/ 1 , . . 1 , . 1 \ 
Q= -W=" (,2 lo ^ + 2m lo S r +2^; 



(22a) 



(22b) 



Proceeding directly with the integration results in the follow- 
ing expression: 



&=( T -1) {cos 0(0^-^4 
-61 



cos 30 



(: 



1 



60r 

3 



16r 7 ^807-V 



# sin 20 



240^^16^ 40r*. 
427 7 3 



__1 
,16r*" 

49 



48r* 
1 

"320r*" 
3 



8r 3T- 12r S 

logr) + 



1 



64r« ' 80r*^ 8T- 2 

_1 3 log r 

80r« 

3 



0+ 



16r 4 

1 . log r> 



1287- 7 " 1 16r) + 



1287- 5 
1 



' ("QGOr 2 
£: sin 40 ( 

Ka ' Mse (i2^-64p-64? 

£?8in 26 (^-327* 

i6? + 240?) +COS 30 ( 48r 48r s ' I61 6 
003 5e (l6T+16^-20>)+^ 20 (n 



512r«^ 
cos 0 
19 5 



_J 3 log r 

"320^ lOr 3 

3 log r\ 
" 64T 2 ^ 



) 



( 



+ 137 



+ 



13 



80r 



3^ 
1 



24? n 
2267 



LA 
Aj 2 ) 



144r». 
23 



+ 



16^2880^ 144r* 



19 , 1 , 11 log r log r log r \ . 
T92r ? " r 40T 5 " r 12^ 2r* t lr r / t 



^sk^f 1 



V 64 48T- 2 288^ 



1 1 1 log r log r\ ■ 

"48r i " t "576r 3 4T- 2 ~4r r /~ r 



^ C03 %2W 



17 

"ger 3 72T- 5 



X s cos 30 



1 3 log r log r log t- N , 
2r 5 ~ t ~ 16r 8r> " t "48r s / /" r 

/— 1 911 , 1 1 log r 5 log r 

V.16r """" " " ~ 



log r logVN r „ . 2g fJE 
ler 6 4r» 8111 20 V576, 

^ cos 0(^-^3)) 



57607 3 " r 1287 J lQOr 7 16r 
1 , 5 log r 



576r 1 36t^ t 32^ 



J28r 1287- 3 , 
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127* 

(23) 



Using the expressions for <fo, <j> u and fa the two-term 
approximation for velocity potential is 



^=U(<h+<t> 1 M 2 +<hM i ) 



(24) 



and from this function the values of velocity at any point in 
the plane may be computed. Of particular interest is the 
evaluation of 

r 00 

for this gives velocity normal to the radius vector of the 
point in question and thus, when r=l, is equal to the velocity 
at the surface of the cylinder. 

Neglecting all terms containing powers of 1/r greater than 
the first, Glauert (reference 1) has given his well known 
result 

y-r Te) = u v 31 * e+ Tr 

and under the same restrictions equation (24) gives 

K KM 2 cos 20 



(-i3>K*H 



2r 4r 
KM* cos 20 , KM* cos 40" 



8r 



16r 



) 



These results are identical to the order of M*. 
Velocity at the surface of the cylinder is 

"2 sin 0 sin 30 



V(lfi) = £/{ 2 sin 0+f +M*[^ 3 2 

2 -^0 + ^]+M{( 7 -l)(f o sin 0- 

11 1 127 

^ sin 30+g sin 50-2 4O X cos 20+ 



|J if cos 40+i|iP sin 0- 



81 



320 



K 2 sin 30- 



Y^K 3 cos 20+^ K* sin 0 



37 

40 



sin 0— 



|j sin 30+^ sin 50- 



887 
"720 



cos 20+ 



H^cos 40+^^8^ 0- 



631 
"960 



K 2 sin 30- 



288 



cos 20+^X* sin 0^)JJ 



(25) 
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Results derived from this equation will be presented in 
terms of pressure coefficient which is defined as 



V-Vo 
1/2*17" 



(26) 



where the zero subscripts refer to free-stream conditions. 
By means of Bernoulli's equation it follows that 



Pm-o — 1 — 



U ' 



(27) 



where P Jf _ 0 denotes the pressure coefficient for incompres- 
sible fluids. If P denotes pressure coefficient for a compres- 
sible fluid obeying the adiabatic law, then 

where M is the Mach number of the free stream and y is 
the ratio of specific heats (1.40 for air). 

As an approximation for P the Glauert-Prandtl result 
(reference 2) is given by equation (2) and the Karman-Tsien 
result (reference 2) is given by equation (1). 



The velocity at the topmost part of the cylinder may be 
found by setting 5=90° in equation (25) and the resultant 
expression is a function of K and M. In figure 1, pressure 
coefficient at this point is plotted against M for K=0. 
As a test for rapidity of convergence the expressions for 
velocity, using only M 1 (one-term approximation) as well 
as M* and M* (two-term approximation), are used. It is 
to be noted that the curves diverge greatly near the critical 
Mach number, but that for smaller values of M, the curves 
derived from equation (25) are together and definitely he 
between the results derived from the Glauert-Prandtl and 
K&rman-Tsien relations. Figures 2 and 3 show the same 
equations applied for i?=l/4 and 1/2, respectively. It thus 
appears from these calculations that the true value of P lies 
somewhere between the approximations applied. On the 
other hand experimental data, as determined from airfoils, 
have shown a better agreement with the Karman-Tsien 
equation than have the theoretical results obtained here for 
the cylinder. 

In figures 4 and 5 the same point on the cylinder is under 
consideration, but circulation is made negative by setting 
K equal to — K and — % in the two cases. As the pressure 
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FiaTOK 1.— Variation of minimum pressure coefficient with Mach number when circulation 

IsO. 

848110— BO 10 



-4.3 



-4.2 



-4J 



t? -4.0 

I 

D 
O 

u 



n i 

i 



-3.9 



-33 



-3.7 



-36 



-35 



T 



— G/auert- Prandff formula 
- Harmon - Ts'ien formula. 



kandt) ; 



One term approximation 

■ Two term approximation 



Circulation - — 
4 




7^ 



-A 



T 



/ 



Critical 

-Mach 

number 



2 .3 
Moch number, M 



Figure 2.— Variation of minimum pressure coefficient with Mach number when circulation 
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Figure 3.— Variation of minimum pressure coefficient with Mach number when circulation 

coefficient gets smaller in absolute value, the theoretical 
data agree more nearly with equation (1). For all the 
calculations, the one-term and two-term approximations 
diverge widely as the Mach number increases. This is to 
be expected for, as has been pointed out by Messrs. G. I. 
Taylor and C. F. Sharman in reference 14, the convergence 
of the series fails when M reaches its critical value. For 
near-critical velocities, several more terms would be re- 
quired to furnish an accurate evaluation of the true potential- 
flow pressure coefficient. 

Figure 6 shows the value of pressure coefficient at all 
points on the surface of the cylinder. These results were 
derived from equation (25) with K set equal to % and at a 
Mach number of %. Crosses on the graph are at positions 
obtained from equation (1) and the circles were determined 
by equation (2). The disagreement at the extreme pressure 
coefficients is again in evidence. 

APPENDIX 

For the integration of the differential equation 
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Fioubk 4— Variation of pressure coefficient with Mach number at topmost point of cylinder 
when circulation b — H rU. 

with boundary conditions 

(!9,.,=°» d ($,..=» < 3 °> 

it is assumed, as in Poggi's method, that the unit circle with 
center at the origin is surrounded externally by a continuous 
distribution of sources such that the ^source strength of an 
element RdRdw is 



j{E, cS)MRdw 



(31) 



Equation (29) may then be interpreted as the equation for 
incompressible flow in such a region. 

The velocity potential of an incompressible fluid at point 
(r, B) due to a unit source at (i2,o) may be calculated by the 
method of images. If this potential is denoted by 0, then 

* = 2^r{! log ^+ R *- 2Rr 008 (*-«)]+ 

^logfp-KR 3 -^ cos (*-«)]+ log |) (32) 
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4>(r,6) =fJJ (log [r>+R?-2Rr cos (8-0,)]+ 
log [^+B 2 -^ cos (0-co)]+ 
2 log ^)/(#,«) BdRdo> (33) 



where the integration extends over the region of the plane 
lying external to the unit circle. 

In the equations under consideration in this report the 
function/ (R, w) is restricted to one of the forms 



sui rrua 
sin mu 



cos ffla 
R' 

COS 771 CO 



logB— gr- logB— ^~ 
m,«^ 1 

As an example of the integration process, take the first case 
listed. Then, set r'= 1/r, which results in 

^,0) = A j" J^log (r«+S a -2iBr cos (0-o>)) ^^BdBdu+ 

J\g(r'*+R?-2Rr' cos{d-<o)ffi^RdRd<>- 



=J l +J a +J 3 



(34) 



Integrating I s first with respect to w shows immediately 
that its value is zero. 
For purposes soon evident T 3 is written in the form 



sin mo. 



RdLRcUo 



(35) 



Since the log R 2 term vanishes, after integration with respect 
to w, the expression for I 3 may be simplified further by the 
substitution 

log [ 1+ (s) " 2 s 008 - 2 SKs)" 003 n( *- w) 

(36) 

since r'<jf?, and 

i r ™ i i 

+sin 710 sin tiw) dRdco — — ^-l — tt sin TntfcZR 

1 1 sin me , 
2mr m m+s-2 

To integrate ii, the region exterior to the unit circle is 
broken into two parts. The first part is a circular ring exter- 
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Fiqubi 6. — Variation of pressure coefficient with Mach number at topmost point of oylindor 
whan circulation is — HrU. 

nal to the unit circle and extending to the-fixed point r, the 
second part is the remaining portion of the'plane and extends 
infinitely. Then 

/l= rJX( log [ 1+ (?J- 2 ? 003 



sm mo) ■ 



R' 



log[l+(£) - 2 (£) cos (e-o^f^RtWd* 

=Ji+J* (38) 
By use of the same series expansion as was previously used, 



-= — sm rnBaR 

m 



= _?mmtf 1 (r »-,+s_ 1) h m-s+2^0 

2mr m m— s+2 v ' 



(39) 
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Fiqube 6.— Pressure coefficient on surface of cylinder for compressible flow with circulation. 



In the same maimer 



r" sin md 
"2m (m+s— 2)r m+, ~ 3 



(40) 



From equations (37), (38), (39), and (40), the solution of 
equation (29), for the case in which 



/(M)=^, 



is consequently 

sin mfl 

* — m(m— 8+2) (m+s— 2) 



_ — sin rofl f 1 1 



360 



4> 



mr- l2m r 2 l0gr ) 



(8—2) ml, , „ , 
jJ=i\ wnen "i+2^s 

when m+2=s (41) 



For the other cases the integration process follows oxactly 
the same procedure. 
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